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ABSTRACT 

We present a new analysis of the relative detectability of dark matter annihilation in the Milky 
Way's eight 'classical' dwarf spheroidal satellite galaxies. Ours is similar to previous analyses in 
that we use Markov-Chain Monte Carlo techniques to fit dark matter halo parameters to empirical 
velocity dispersion profiles via the spherical Jeans equation, but more general in the sense that wc do 
not adopt priors derived from cosmological simulations. We show that even without strong constraints 
on the shapes of dSph dark matter density profiles (we require only that the inner profile satisfies 
— limr_>.o dlnp/dlnr < 1), we obtain a robust and accurate constraint on the astrophysical component 
of a prospective dark matter annihilation signal, provided that the integration angle is approximately 
twice the projected half-light radius of the dSph divided by distance to the observer, a^^t ^ '^fh/d. 
Using this integration angle, which represents a compromise between maximizing prospective flux and 
minimizing uncertainty in the dSph's dark matter distribution, we calculate the relative detectability 
of the classical dSphs by ground- and space-based 7-ray observatories. 

Subject headings: galaxies: dSph, dark matter profile — dark matter: indirect detection — methods: 
MCMC 



1. INTRODUCTION 
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tents (Ly - lO^-^Lyo), the Milky Way's dwarf 
spheroidal (dSph) satellite galaxies have become tar- 
gets of interest in searches for high-energy ph otons 
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as-yet unidentified particle(s) that constitute cosmo- 
logical dark matter. Such an 'indirect' detection 
of dark matter has thus far remai n ed elusive (e.g. , 
Abdo et all 120101: lAcciari et all [20T0I: lA haronia n et al 
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201 in . but recently-commissioned and planned instru- 



ments w ill soon explore the sky with unprecedented sen- 
sitivity (jAtwood et al.l[2009l : ICTA Consortiumll2010f ). 

Any prospective signal depends on the nature of dark 
matter and its distribution in the emitting object. Stellar 
kinematic data can provide information about the distri- 
bution of dark matter in dSphs that, combined with some 
future observation of annihilation/decay products, would 
then amount to information about the nature of dark 
matter. High-quality stellar kinematic data sets now ex- 
ist for each of the Milky Way's eight 'classical' dSph 
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have used these data to predict annihilation signals 
(e.g., iBringmann eTall [20091 : iPieri et al.l [200l) Sev- 
eral published calculations fe.g.. iStrigari eraU [20081 
iMartinez et al.| [2009: Kuhlcn 201^ of indirect detection 
signals from dSphs begin by adopting 'cosmological pri- 
ors', effectively assuming that dSphs occupy dark matter 
halos that have structural parameters (e.g. virial mass, 
scale radius and asymptotic behavior of the inner density 
profile) characteristic of halos produced in dissipation- 
less N-body simulations o f galaxy form ation that assume 
dark matter is 'cold' f e.g..lNavarro. Frenk & White 19971: 
Navarro et all [200l iKuhlen et al.l 120081: ISpringel et al.l 
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Given a lack of direct observational constraints on the 
density profiles of dSph dark matter halos, here we pro- 
vide a more general analysis with the goal of examining 
the robustness and veracity with which a Jeans analy- 
sis of dSph stellar-kinematic data can constrain the as- 
trophysical component of a dark matter annihilation sig- 
nal. We adopt minimal assumptions about the structural 
properties of dSph dark matter halos and we test our 
method on artificial data drawn from physical distribu- 
tion functions. In this context we identify the optimal in- 
tegration angle that balances interests in maximizing flux 
and minimizing uncertainty in the astrophysical compo- 
nent of a prospective annihilation signal. Finally, we 
compare the classical dSphs in terms of the relative de- 
tectability (absolute detectability depends on unknown 
particle physics) of their dark matter annihilation signa- 
tures. In subsequent contributions we revisit this topic 
in greater detail, considering the impact of (sub-) sub- 
structure and providing a thorough discussion of the in- 
strumental response (Charbonnier et al. in preparation). 

2. J-FACTOR AND CONFIDENCE INTERVALS 
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The differential 7-ray flux (units cm~^ s^^ GcV^^) 
received on Earth in solid angle Ail from dark matter 
annihilations is given by 
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where is the dark matter particle mass, (av) is the 
velocity-averaged annihilation cross section, dN^/dE^ is 
the energy spectrum of annihilation products and 
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This 'J- factor' represents the astrophysical contribution 
to the signal and is specified by the integral of the 
squared dark matter density, O), over line of sight / 
and solid angle fl. 

In order to allow for a broad range of dark matter 
density profiles, we adopt the |Zhao ( 1 996h genera liza- 
tion (i.e., the a, /g, 7 mod e l) of t he iHernquistI (jlOOO") and 
iNavarro. Frenk fc Whiti (|1996l 'NFW hereafter) pro- 
files, 



P{r) = Ps{— 



-7 r 



1+ - 



(3) 



We do not consider subclumps in the dSphs, as these 
substructures ar e not expected to con tribute significantly 
to the J- factor (jSpringel et al.|[2008l Charbonnier et al., 
in preparation). 

2.1. Toy-model behavior 

The J-factor depends primarily on the inner logarith- 
mic slope 7 = — lim^-^o dlnp/dlnr, scale radius r^, and 
normalization ps of the dark matter halo. Only two of 
these three quantities are independent for the classical 
dSphs, such that masses within 300 pc c an be reason- 
ably w ell-constrained at M300 ~ IO^Mq (jStrigari et al.l 
I2008al) . Here we find that J-factors for the classical 
dSphs depend only weakly on the exact profile as long 
as the inner profile has slope 0<7<1(7 = cor- 
responds to a 'core' of constant density, whereas 7 = 1 
corresponds to an NFW 'cusp'). This independence can 
be understood using the point-like approximation for- 
mula for J, expressed as a function of dSph distance d 
and integration anglcQ aint (with AQ = 2n{l — cosaint)) 
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The shape of the integrand in Equation 2] (using a log- 
arithmic integration step) is plotted in Fig. [T] for several 
(a,/3,7) = (1,3,7) dark matter profiles (i.e., NFW with 
a free inner slope). Each is normalized to the value at 
r = rs/2, a radius that is expected to be of the same 
order of magnitude as the typical radii of stellar orbits. 
For r > rsi the contribution of the integrand quickly van- 
ishes in all cases as the integrand drops at least as fast as 
r""*. For 7 = 1.5, the integral diverges in the inner parts 
unless one imposes a cutoff (such a cutoff is physically 

® To avoid confusion with halo parameter a appearing in Equa- 
tion [3] wc use subscripts to denote integration angles, e.g., ajnt 
and Oc- 
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Fig. 1. — Integrand appearing in Eq. Q to calculate the J-factor 
for (1,3,7) dark matter profiles. For comparison we show the in- 
tegrand calculated for the Einasto profile. In this representation, 
■^point-liko is directly proportional to the area under each curve. 



motivated by the need to balance rates of an nihilation 
and gravitational infall. iBerezinskv et al.|[T992D . 

Otherwise, all the integrands of profiles with < 7 < 1 
arc peaked and have similar areas under their curves. For 
1 < 7 < 1.4, the integral is larger by a factor of a few, 
with the exact value depending on the choice of normal- 
ization. Figure [U also show s the integrand calcula ted for 
the 'Einasto profile' (e.g., lEinasto fc Haudlll989( ) that, 
with respect to NFW and by virtue of its variable inner 
slope, provides slightly better fits t o simulated cold dark 
matter halos (jNavarro et al.l |2004[) . The three orange 
curves in Figure [T] corres pond to the Einasto p rofiles (in- 
dices n = 4, 6, 8) used bv INavarro et al.l (|2004[ ) to fit the 
density profiles of low-mass cold dark matter halos. The 
integrands from Equation 2] corresponding to the plotted 
Einasto profiles (chosen here to have — dlnp/dlnr — 2 
at the scale radius of the plotted NFW profile) behave 
similarly in Figure [1] to the a, /?, 7 models plotted with 
7 < 1.5. 

2.2. Jeans/MCMC analysis of classical dSphs 

Nearly all dSph stars with measured velocities lie 
within a few half-light radiQ of the dSph center (r^ ^ a 
few xlO^ pc for the Milky Way's classical dSphs). Var- 
ious published analyses have used the Jeans equation to 
show that th e observed flatness of dSph velocity disper- 
sion profiles (jWalker et al.ll2007| ) leads to a constraint on 
M{rh) — the mass enclosed within a sphere of radius — 
that is insensitive to assumptions about cither anisotropy 
or the structural parame ters of the dark ma t ter halo as 
speci fied, e.g., by a, /3, -7 (I Walker et"all l2009l: IWolf et al.l 
l2010l : lAmorisco fc Evansll2011|) . For an indirect detection 
experiment, locating the optimal integration angle, de- 
noted ac, on a resolved target amounts to finding a com- 
promise between maximization of signal (in the presence 
of background) and minimization of model-dependence 
in the corresponding astrophysical factor Jq,^ — and hence 
minimization of the uncertainty on the derived values of 
and {crv). For the Plummer models often projected 
to fit dSph surface brightness profiles, 80% of the pro- 



For consistency with IWalker et all I I2009I) we define as the 
radius of the circle enclosing half of the dSph stellar light as seen 
in projection. Elsewhere this radius is commonly referred to as the 
'elfective radius'. 
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Fig. 2. — J median value and 95% CLs from four independent 
MCMC analyses of Carina, holding the parameter 7 fixed at the 
values shown. The arrow identifies ajnt = etc ~ Sr^/d. 



jected stellar mass lies within 2rft. An integration angle 
aint = etc ^ 2rh / d therefore represents a reasonable com- 
promise between flux collection and robust estimation of 
the J-factor. 

We verify the robustness of Jeans estimates of J^^ — 
i.e., the J-factor corresponding to integration angle 
Q^int = otc ^ 2r/,/(i — by repeatin g the Markov - Chain 
Monte Carlo (MCMC) analysis of I Walker et"all (pOOl 
WOO hereafter) for the Milky Way's eight classical dSphs. 
This analysis assumes virial equilibrium, spherical sym- 
metry, constant velocity anisotropy, a dark matter den- 
sity profile given by Eq. ([3]), and a stellar luminos- 
ity profile given by a Plummer model, L{r)/Ltot = 
(r/r;i)^(l + [r/r/i]^)^^/^, with half-light radius derived 
from the s tructural parameters originally published by 
llrwin fc Ha tzidimitriou (1991). We use the same velocity 
disper sion profiles displayed in Figure 1 of I Walker etall 
()2009D (estimated under the assumption that the veloc- 
ity distribution within each bin is Gaussian, thereby dis- 
carding dynamical information that might be containe d 
in non-Gaussian distributions; e.g., iLokas et al.l l2005h . 
but our current analysis differs in the following ways. 
First, we allow the outer logarithmic slope of the dark 
matter density profile to remain a free parameter, /3, for 
which we take a flat prior over the range 3 < ^ < 7 (this 
change has no significant impact on our results) . Second, 
while we eventually allow 7 to vary freely between and 
1 (section [3]), we perform several initial analyses of each 
dSph in which we fix the inner slope of the dark matter 
density profile at specific values of 7 = (a constant den- 
sity core), 7 = 0.5 (a mild cusp), 7=1 (NFW cusp) and 
7 = 1.5 (a ste ep cusp). Th i rd, aft er the MCMC chains 
converge (see I Walker et al.l (|2Q09( ) for complete details 
of the MCMC procedure), wc perform a numerical inte- 
gration of Eq. ([2]) that yields a value of the J-factor for 
each set of parameters present in the chains. The distri- 
bution of these values then represents the observational 
constraint on the J-factor, given the modelling assump- 
tions. 

Figure [2] displays the median J-factor and 68% con- 
fidence intervals (CIs) from our MCMC chains for the 
Carina dSph, where we have used integration angles 
Qfint between 0.01 and 5° (a, /3, r^, and the veloc- 
ity anisotropy parameter free) . In order to make explicit 
the effect of strong assumptions about the inner slope 



of the dark matter density profile, we show results from 
our four initial MCMC analyses with the slope held fixed 
at values 7 — 0.0,0.5,1.0,1.5. Small integration angles 
Q^int Qfc '2rii/d in Fig. [2] correspond to r <C ?'s in 
Fig.lll so that the different values for the adopted 7 result 
in different median values and CIs for the J-factor. How- 
ever, for integration angles a-mt = cic Sr/j/d, the me- 
dian J-factors and corresponding CIs converge to similar 
values regardless of whether we adopt 7 = 0. 0,0. 5 or 1.0. 
Adoption of a larger inner slope (e.g. 7 = 1.5 in green) 
gives a larger J-factor even as evaluated at aint = etc, as 
expected based on the point-like approximation shown in 
Fig.[TJ The other dSphs (not shown) exhibit the same be- 
haviors, although with narrower/broader confidence in- 
tervals for dSphs in which more/less kinematic data are 
available. 

2.3. MCMC analysis of artificial data 

In order to examine its reliability, we apply our 
Jeans/MCMC analysis to artificial sets of stellar posi- 
tions and velocities drawn from distribution functions 
of the form L~^^f{e). Here e is energy, L is angular 
momentum and b is the velocity anisotropy parameter 
b = 1 — (Tj /fT^ , where and are second moments 
of the velocity distribution in radial and (one of the) 
tangential directions, respectively. For 6 0, these 
models give stellar distributions wi th constant velocity 
aniso tropy (for a discussion see, e.g ., lBinnev fc Tremaing 
120081 ). Following [Cuddeford' (1991|), once the halo model 
and stellar density are specified, an Abel inversion is used 
to determine the distribution function numerically. We 
generate a grid of halo models with 7 = 0.1,0.5,1.0; 
Th/rs = 0.1,0.5,1.0 and (3 = 3.1. For each model we 
consider isotropic (b ~ 0) and anisotropic distributions 
with b = 0.25 or b = —0.75. We also generate a grid of 
models with a = 1.5, /3 = 4.0. All test halos are nor- 
malized to have mass ~ IO^Mq within 300pc. Stellar 
positions and velocities are generated directly from the 
full stellar distribution function using an accept-reject 
algorithm, and subsequently all velocities are scattered 
randomly to mimic observational errors of ±2 km s~^. 
All artificial samples have 1000 stars, similar to the typ- 
ical data set available for the classical dSphs. 

In our analysis of the artificial data sets, as in our 
final analysis of the real dSph data (section [Sj , we al- 
low 7 to vary freely between < 7 < 1 (effectively 
adopting a flat prior over this range) , reflecting the facts 
that 1) published kinematic analyses d o not place strong 
constraints on this parameter (e.g. iKoch et al.l 120071 : 
IWalker et al.l 120091) . and 2) current cosmological simu- 
lations suggest that 7^1 f or the sub halos most readily 
associated with dSphs (e.g.. lSpringel e t al. 200^. For 72 
artificial dSphs, Fig. [3] plots JMCMc/t^truc, the ratio of 
the measured (median and 95% CIs) to the known input 
value of J, as a function of aint/ac- For most of the test 
models, our estimates of Jmcmc (o^c) he within a factor 
of two of the correct values, even for models with ex- 
treme anisotropy. However, for the most steeply cusped 
(7 = 1.5) test models that violate our assumption that 
the inner slope satisfies < 7 < 1, we under-estimate 
the J-factor typically by a factor of > 10 (green curves 
in Figure While the adoption of a less restrictive 
prior (e.g., allowing < 7 < 2) would avoid this po- 
tential bias, such an improvement would come at the 
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Fig. 3. — From artificial data sets, ratio of the measured J-factor 
to the true J-factor, as a function of ai^i/ac- Lines are color-coded 
with respect to the true value of the slope of the inner density 
profile of the input model. For these analyses, as for the real dSph 
data (section [Sj, we adopt a uniform prior over < 7 < 1 for 
the inner slope. Solid lines indicate median values of the J-factors 
obtained in our MCMC chains, while dashed/dotted lines indicate 
lower/upper bounds on 95% CIs. 



cost of larger uncertainties for estimates of the J-factor 
in all halos. Since c urrent cosmological simulations (e.g., 
iSpringel et al.ll2008f) and indirect arg uments based on ob- 
servations fe.g.. IKlevna et al.ll200a IGoerdt et "all 120061: 
IGilmore et al.l 120071 ) suggest that dSph dark matter ha- 
los have 7 ^ 1, we report results based on our more 
restrictive assumption that < 7 < 1. Figure [3] illus- 
trates that our estimates of the J-factor are as reliable 
as this assumption. 

3. DETECTION PROSPECTS 

Table [1] and Figure S] indicate estimates of Jq^ we ob- 
tain for the real dSphs from our Jeans/MCMC analysis 
while letting 7 vary freely between < 7 < 1 (i.e., the 
same analysis tested using artificial data). The classi- 
cal dSphs all lie well above the Galactic plane and hence 
are prone to similar levels of diffuse 7-ray background. 
Adopting the background expected for each dSph from 
the standard Fermi diffuse model of the Milky Way (S. 
Funk, private communication), we find < 7% rms vari- 
ation from dSph to dSph in estimated background flux. 
We neglect these small background variations in the dis- 
cussion that follows. 

For a signal-limited detection (plausible in the case of 
the Fermi-LAT) Jq,^ provides an appropriate figure of 
merit for selecting dSphs as dark matter annihilation 
targets (fill ed circles in F i g. [H). In that case we find, 
similarly to iStrigari et al.l (|2007[ ). that Draco and Ursa 
Minor are among the best targets in terms of having the 
largest median-likeli hood estimates for Jp „. Somewhat 
at odds with that of IStrigari et al.l (j2007l ). our analysis 
places Sculptor on approximately equal footing with Dra 
and UMi and perhaps even above UMi given the large 
error bar associated with that galaxy (of the three, UMi 
has the smallest available kinematic data set). This dis- 
crepanc y likely follows from different choices of prior for 
p{r), as' lStrigari et al.l (|2007j ) assume NFW profiles. 

On the other hand, in the presence of background, ob- 
jects with ac significantly larger than the instrumental 
point-spread-function are disfavored. Taking 0.1° as an 
optimistic estimate of the angular radius containing 80% 
of the PSF (close to the optimum integration radius for 
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Fig. 4. — Astrophysical factor as a function of the distance 
to the dSph. Filled circles, solid and dashed error bars corre- 
spond respectively to the median value, 68%, and 95% CIs on 
Ja^ (where ac is given in Table [TJ. Empty circles correspond to 

= J.JVi + K/o.i°)2. 



a point-like source in the background-limited case) of fu- 
ture ground-based 7-ray detectors in the approp riate en- 
ergy range (0(mlc^/5). \CTA Consortiumll2010l) . we use 

~ Jq^/-\/1 + (q!c/0.1°)^ as a figure of merit in the 
background- limited cascQ (empty circles in Fig. |4]). A 
dominant background significantly changes the ranking 
of the dSphs, with Leo II now among the best targets. 
Note, however, that Leo II has the largest error bar as it 
also has the smallest available kinematic data set of the 
eight galaxies considered here. 

The J-factors estimated here do not reach the val- 
ues required for dark matter detection with existing and 
near-future instruments in the most conventional parti- 
cle physics scenarios. Moreover, most are smaller than 
the values predicted for the Galactic Center (GC). To 
illust rate the last point, we adopt an Einasto profile 
(e.g., iNavarro et"aIll2004D for the Milky Way's smooth 
dark matter component. We then estimate Jgc ~ 
2 X 10"M|kpc-^ toward the GC for ai„t = 0.01°, 

0.1°. How- 



and Jgc ~ 1 x lO^'^ Mg kpc ' for ai„ 
ever, the presence of st rong astrophysical back grounds 
on scales of arcminutes ("Aharon ian et al.l |2004[) , tens of 
arcminutes (lAharonian et al., .20(]6() and tens of degrees 
(|Su et al.ll2nmr will obfuscate any genuine dark matter 
signal coming from the direction of the GC. In contrast, 
the old, gas-free stellar populations of dSphs provide few 
sources of point-like or diffuse 7-ray emission. On these 
grounds dSphs provide a favorable alternative to the GC 
as targets for indirect dark matter searches using instru- 
ments with resolution a-mt ^ 0.1°. 

4. CONCLUSION 

We have shown that our Jeans/MCMC analysis, in 
which we have allowed the inner slope 7 to remain a free 
parameter with possible values between < 7 < 1, esti- 
mates the astrophysical J-factor to within a systematic 
uncertainty of a factor of < 3 for the Milky Way's classi- 
cal dSphs, so long as 1) the integration angle is chosen to 

* We obtain this figure of merit by convolving the dark mat- 
ter signal with the PSF while assuming that both are Gaussian. 
Charbonnier et al. (in prep.) study the effects of angular reso- 
lution and demonstrate that within the 80% containment radius, 
this approximation is accurate to within a few percent. 
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TABLE 1 

Classical dSphs sorted according to their distance I|MateciII1998I '). The other columns correspond to twice the half-light 
radius, the optimum inteciration ancile, the median log]^g(j) with 68% (95%) confidence intervals for that ancile, and 

background-corrected logjQ(J^'2) (see text). 
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be Q!int = ac ^ '2rh/d and 2) the actual dark matter halo 
has 7 < 1. For cuspier (7 > 1) profiles, Figures [T] and 
[3] indicate that our analysis under-estimates the J- factor 
systematically (e.g., by a factor of > 10 if the input halo 
has 7 > 1-5), but such steeply cusped profiles are neither 
supported by observations nor motivated by current cos- 
mological simulations. For the family of Einasto profiles 
applie d to dSph-like dark matter halos by iNavarro et alj 
dloP), our method over-estimates J by a factor of < 2 
(Figure [T]). We conclude that our analysis provides esti- 
mates of the J-factor that are reliable for the dark mat- 
ter halo models (e.g.. iMerritt et ani2006L and references 
therein) most relevant to dSphs. 

We consider Figure |3] to provide a good starting point 
for ranking the classical dSphs as targets for indirect de- 
tection of dark matter via annihilation. For a given parti- 
cle physics scenario the highest median-likelihood values 
of Ja^ correspond to the largest expected fluxes. In the 
background-limited case, appropriate for searches with 
future ground-based 7-ray observatories such as CTA, 
the highest values of J^^ correspond to the largest ex- 
pected statistical significance. Figure |4] illustrates the 



large differences in the fractional uncertainty in J^^ be- 
tween different dSphs. This uncertainty also forms a cri- 
terion for target-selection. In the end, those dSphs with 
simultaneously the largest and most tightly constrained 
J-factors have the greatest potential to provide useful 
constraints on the particle physics of dark matter. 
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